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Abstract 
Introduction. The problem of critical loads of a compressed orthotropic rectangular plate on an elastic base was 


considered. The following orthotropy parameters were set for the plate: Poisson coefficients, Young's modules for the 
main directions, and the shear modulus of the plate material. The components of the compressive load were uniformly 
distributed along two opposite edges of the plate and acted parallel to the coordinate axes. The edges of the plate were 
loosely pinched or pivotally supported. Cases were also considered when two parallel edges of the plate were free from 
loads, and the other two were freely pinched or pivotally supported. 

Materials and Methods. The problem was studied on the basis of a system of nonlinear Karman-type equilibrium 
equations. The critical values of the load parameter were determined from a linearized problem based on a trivial 
solution. At the same time, the variational method in combination with the finite difference method was used to solve 
the boundary eigenvalue problem. 

Results. The problem was reduced to solving a parametric linear boundary eigenvalue problem. In case of boundary 
conditions of a movable hinge support, exact formulas of eigenvalues and eigenfunctions were given. While in case of 
free edge pinching, a variational method was used in combination with a finite-difference method, and a computer 
program for solving the problem was built. It was established that one or two eigenfunctions expressing the deflection 
of the plate could correspond to the critical value of the compressive load parameter at which the stability of the 
compressed plate was lost. The results of numerical calculations of the critical values of the compressive load at 
different values of the orthotropy parameters were presented, and graphs of the corresponding equilibrium forms were 
constructed. For the case of a long orthotropic plate on an elastic base, it was established that the main term of the 
asymptotic expansion of the solution to the linear eigenvalue problem was determined from the problem of critical loads 
of a compressed beam on an elastic base with an elastic modulus that coincides with the elastic modulus of the plate in 
the longitudinal direction. 

Discussion and Conclusions. The problem of critical loads of an orthotropic plate compressed in two directions lying 
on an elastic base was investigated. As the compressive load component increased along one direction, the critical value 
of the load compressing the plate along the other direction decreased. If an orthotropic plate was compressed by a load 
along a direction that corresponded to a greater bending stiffness, then the critical value of the loss of stability was 
greater than the critical value of the compressive load acting along the direction of a lesser bending stiffness. The 


presence of an elastic foundation increased the bearing capacity of the compressed plate. 
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Introduction. L. M. Zubov in [1] derives the equilibrium equations of an elastic plate containing sources of internal 
stresses in the form of dislocations and disclinations, which are a modification of the Karman equations. In this work, 
the problem of bending a thin plate (membrane) under the action of internal stresses caused by defects is also solved. 
In [2-5], the theory of dislocations and disclinations is used in the study of the behavior of plates and shells. In 
article [6], a general solution of the oscillation equation of a rectangular orthotropic plate with free edges is constructed. 
In [7, 8], asymptotic solutions of problems of statics and dynamics of narrow plates are constructed. 

In [9], the problem of equilibria branching under the influence of a small normal pressure of a rectangular plate 
compressed in one direction with dislocations and disclinations is considered. The case is considered when the 
longitudinal edges of the plate are free from loads, and the other two edges are fixed. In [10], an asymptotic solution is 
provided for the problem of critical loads of a compressed narrow plate using a small parameter expansion of the 
relative width of the plate. In [11], the branching problem of the equilibria of an elastic rectangular plate with internal 
stress sources is investigated for the case when compressive loads are unevenly distributed along the edges. It is 
established that with even forms of the incompatibility function and even forms of distribution along the edges of 
compressive forces, the presence of a small normal load does not reduce the bearing capacity of the plate. 

In work of Morozov N. F., Belyaev A. K., Tovstik P. E., Tovstik T. P. [12], an asymptotic derivation of two- 
dimensional equilibrium equations of a thin elastic inhomogeneous plate made of a general anisotropic material is 
presented. Conclusions about the forms of loss of stability of a compressed plate on an elastic soft base are formulated 
in [13]. In [14], the bending stiffness of a thin elastic multilayer plate with transversely isotropic layers is investigated. 
Two models of accounting for the effect of transverse shear are compared. 

In [15], a model of thermoelastic viscoplastic deformation of a composite that is cross-reinforced with continuous 
fibers in arbitrary directions is constructed. 

In [16], within the framework of a direct approach to plate theory, the problem of a hyperelastic plate with 
inhomogeneously distributed initial stresses is considered. The plate is considered as a material surface with five 
degrees of freedom (three displacements and two rotations). In [17], the stability of a uniformly compressed circular 
two-layer plate with an initially compressed or stretched layer is analyzed within the framework of nonlinear elasticity. 
A model of incompressible neo-Hookean material is used for determining the material ratio. The analysis of the 
dependence of the resulting critical stresses on the initial deformations and stiffness parameters is presented. In [18], the 
bending of a three-layer plate with surface and interfacial stresses is investigated, the theory of plates with first-order 
shear deformations and the Gurtin-Murdoch model of surface stresses are considered. The dependence of the plate 
stiffness parameters on the surface elasticity modules is analyzed. In [19], the problem of instability of a three-layer 
nonlinear elastic rectangular plate with a prestressed middle layer is considered. 

In this paper, the influence of the parameters of the plate material and the elastic base parameter on the critical loads 


of the loss of stability of a compressed elastic rectangular orthotropic plate on an elastic base is investigated. 
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Materials and Methods. Let us consider an orthotropic rectangular plate lying on a linearly elastic base, which is 


compressed by forces P and Q along the X and Y axes. The equilibrium equations can be written as [7]: 


D,0,. W + 2D,0;,.0, W + D,0,W+KW =[W,F], 


2 (1) 
 otpy| £—-™ leeatrs + atr=—[w,w]. 
E, G E£, E, 2 
Let us write the boundary conditions in the form: 
W=0,W=0, 6,0,F=0, 60;F=-Pat |X|=a/2, bs 
W=0,W=0, 6,0,F=0, 0,F=-Qat |¥|=b/2. 
W=6,W=0, 6,0,F=0, 0,F=-Pat |X|=a/2, és 
W=0,W=0, 0,0,F=0, 0:F=-Qat |Y|=b/2. 
W=0,W=0, 0,0,F=0, O0;F=—P at |X|=a/2, 4) 
0, W +v,0;,W = 0, W +(2-v,)0,0,W =0,0,F = 6, F =0 at |Y|=b/2. 
W=6,W=0, 0,0,F=0, O0;F=—-P at |X|=a/2, 6) 
0, W +v,0;,W = 0, W +(2-v,)0,0,W =0,0,F = 6,F =0 at |Y|=b/2. 


[W,F] = 05, W0,F+ 0, W0;,F-26,0,W0,0,F, D,=Eh’/(12(1-v,v,)), i=12; D, =Gh’/12, D,=D,v,+2D,, 
E,v, =E,y, . 


Here v,, v,, E,, E,, G— Poisson's coefficients, Young's modules for the main directions and the shear modulus 
of the plate material; D,, D, — bending stiffness in the main directions; D, — torsion stiffness; the function 
W(X, Y) expresses deflection, and F(X, Y) denotes the stress function; 4 — plate thickness, KW — base reaction; 


K — modulus of subgrade reaction; the origin of coordinates X, Y is in the center of the plate and the axes are parallel to 
its edges; 0 — partial derivative of the n-th order with respect to the variable X. We assume that the edges of the plate 


are loosely pinched (2) or pivotally supported (3). In this case, the compressive forces P and Q are evenly distributed 


along the edges X =+a/2 and Y=+b/2 respectively. The cases are also considered when two parallel edges of the 
plate Y =+b/2 are free from loads, and the other two are freely pinched (4) or pivotally supported (5). In these cases, 
the compressive load is applied only to the edges X =+a/2. 


Y* x 
Let us assume that F=®-—P 5 Q a Then, for the function ® , homogeneous boundary conditions are valid 


®,, =D, =0 at [X|=a/2, O,, =®,, =0 at /¥|=b/2. These conditions can be replaced by equivalent boundary 
conditions ®=®, =0 at |X|=a/2; ®=@,=0 at |Y|=b/2. Let us move on to dimensionless variables by 


formulas: 


KD, poling 


. _qD, . 
Ae ae 


b 


X=ax; Y=by; 6 » . @(X,Y)=Df (xy); W(X,Y)=w(xy)h; K= 
a 


Then system of equations (1) and boundary conditions (2), (3) can be rewritten as 
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ees +2¢,8°O?0?w + c,0'w + 8‘kw = 8? [w,f]—p3°d2w —qd?w, é 
S*d%f +2c,02025"f + c,04f = -ac,8° [w, w]. 
[w.d,w.f,d,f]},_,, =05 [w, Aw fd, |, = (7) 
[w.dw.f.0f |, =05 [ w.2w. £0, | Lag (8) 
[w.d,w.f,d,f],_,, =95 [O,w +v,0,w,d,w +(2—v,)0,0,w, £,8,f ] =O, (9) 
[w.aw.f.aF |, =0; [@w+v,dw, dw +(2-v, )a2d,w, a = (10) 


For any values of the parameters k, p, g the nonlinear boundary value problem (6) with any of the conditions (7)— 


(10) satisfies a trivial solution (w..f,) = (0,0). 
Let us denote by E’ the set of vector functions f = (f.6,) with a finite norm, which is determined using the scalar 


product: 


(f.8) =[J(fieitfs.)dxdy; f =(f,,f,).¢=(g,.82)- (11) 
Q 
By E! we denote the space of infinitely differentiable vector functions u=(w,F), v =(w,,F) in the domain 


Q= {(x,y):|x| e1/2,, 


domain Q . We define the norm in E' using the scalar product 


(u,v). = | ie Sa (12) 
ivjs4 2 


y| <1/2} , which satisfy one of the boundary conditions (5) or (6) on the boundary OQ of the 


Ox'dy! , Ox'dy! 
The boundary value problem (6) with one of the boundary conditions (7)—(10) can be considered as a nonlinear 


operator equation 


M,u=Ih, u=(w,f)eE', (13) 
where 
5°O¢w + 2¢,8°0.0,w +¢,0)w + kw + p60. w + g0,,w a 3°[w,f] as 
— 5 U= : 
*  8tatf +. 2¢,0°628"F +0,0°F —ac,5°[w,w] 


Here M, — linear, and II — nonlinear operators, since 
[w.f] — 0. w0,f +0,w0;f —20,0,w0,0,f . 

From the results of the works of I. I. Vorovich and N. F. Morozov, it follows [9] that the operators M, and II act 
from the space E! to E’. Let k = 0, p = 0, q = 0. Then equation (9) has a unique trivial solution u, = (w.,f.) = (0,0). 
With small changes in the values of the parameters p, q and k the trivial solution remains the only one. If, for any 
small values of a number ¢>0 and the fixed value k there are parameter values p, q, that satisfy the inequality 
[P—Po|+|a—q|<« and for which there are at least two solutions to equation (9), then the point (py, qy,) is called the 


branching point (or bifurcation point) [9] of this equation. It is shown in [9] that in order to find the bifurcation point 


(Po. qo) , it is necessary to solve the eigenvalue boundary value problem, which is a linearized equation on a trivial 


solution 


Mu =0. (15) 
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The eigenvalues (Po. qo) are the critical values of the parameters of the compressive load (p. q) . Considering (13) 


and (14), equation (15) can be written as the following eigenvalue boundary value problem: to find values of the 
parameter p, such that for fixed values of the base parameter (modulus of subgrade reaction) k and the load parameter q 


along the y axis, the following problem has a nontrivial solution: 


3°Olw + 2c,8°0.0,,w + c,0,w + kw + q0,w =—ps’d-w , (16) 
[w.d,w]}, y= 9; [w.3,w =0, (17) 

[ w.a2w | pein = 83 [ w.2w | ano (18) 

[w.2,W],. = 0% [Owtvd.w,d,w+(2—-v,)OV,w, f,8,f ]) =0, (19) 
[ w. aw | wa 73 [ ow +v,d,w, dw +(2-v, )d20,w, ,8,f ] =0. (20) 


In this case, the critical value of the loss of stability of the compressed plate corresponds to the minimum eigenvalue 
of the parameter p, and the eigenfunction corresponding to this value shows the form of equilibrium after the loss of 


stability. For the case of the hinge support of the edges (18), exact solutions can be written for equation (16): 


(5‘m* +2c,m’n* +o,n*) x! —qn’n? +k ine 1 
: > Wan =Ssin7m ee sin ys |: (21) 


Pin nm 
where m and n are arbitrary natural numbers. 
From the analysis of formulas (21), it follows that with different ratios of the plate edge lengths 6, the physical 


parameters of the orthotropy of the plate material c,, c,, the parameter g of the compressive load along the y axis, and 


the modulus of subgrade reaction k (characterizing the reaction of the base), one or more half-waves in the form of 


equilibrium (21) can correspond to the critical load value p, = MIN Pian In this case, the so-called buckling by two 
proper forms is possible, when, for example, p, = MIN Py =P; =P2,- Such cases are listed below. To solve the 


problem (16) with boundary conditions of free pinching of the edges (17), a variational method is used in combination 
with a finite-difference method, which was justified in [9, 10] and is easily transferred to the case of an orthotropic 
plate. At the same time, the computer program developed for this task was tested on the solution of the 
problem (16), (18). In addition, a computer program has been developed for solving the variational method in 
combination with the finite difference method of the problem (16) with each of the boundary conditions (19), (20). For 
these cases, test examples of the corresponding boundary value problems have been constructed. 

In [7], in the problem of asymptotic integration of the equation of small oscillations of a long (b << a) rectangular 


orthotropic plate, an expansion of the solution by small parameter degrees 6=b/a is constructed and it is shown that 
the main expansion is determined from the equation of vibrations of a beam with an elastic modulus E=E,. Repeating 


the reasoning of works [7] and [10], we will look for a solution to problem (16) with boundary conditions (19) or (20) in 


the form of series: 


p=) p,8", w= ow," . (22) 
i=0 i=0 
It can be shown that the main term of the asymptotic expansion (22) of the solution to problem (16) with boundary 


conditions (19) or (20) is determined from the problem in dimensionless variables about the critical loads of a 


compressed elastic beam on an elastic base with freely pinched or movably hinged edges: 


(1-v,) wi'+ kw, +powy =0, (23) 


Peshkhoev I. M., et al. Critical Loads of Uniformly Compressed Orthotropic Rectangular Plate on an Elastic Base 


WwW, =Wo =0,~ (24) 

WwW, =Wo =0. (25) 

Here the prime symbol indicates the ordinary derivative of a function depending on the variable x. It is possible to 

construct an algorithm for determining all coefficients of the terms of the expansions (22), but the authors will not do 
this in this paper. 

When moving to dimensional variables from (23)—(25), we obtain the problem of critical loads of a compressed 


beam with an elastic modulus E=E,. The eigenvalue problem (23) with boundary conditions (25) has an exact solution 


Pp, =2'm + 


a ; w,, (x) =sin(zmx). (26) 

In the case of boundary conditions (24), problem (23) is solved by the variational-difference method. 

Research Results. Numerical calculations of the critical values of the load parameter p compressing along the 
direction of the x axis with fixed values of the modulus of subgrade reaction k and fixed values of the load parameter q 
compressing along the direction of the y axis and the specified values of the orthotropy parameters of the plate material 
are carried out. The plate equilibrium forms corresponding to critical loads are constructed. 

Example 1. Let us consider an example of a square orthotropic plate from work [16] "plywood of the first type" — 


three-ply: if you direct the x axis along the fibers of the face (upper and lower of the three layers), then Young's 


modules along the fibers and across, respectively, are equal to E, =1.4-10°kg/cm* and E, = 10" kg/cm’; Poisson 


coefficients — v, =0.46 and v, =0.46/12; shear modulus for the main directions of elasticity G=1.2-10°kg/cm’. 


b D, £E 1 

Then, for the problem in dimensionless variables (16), we obtain 6 1s 2 2 0.083 ; 
a D, £E, 12 

c, =— =0.207 taking into account the above formulas for the transition to dimensionless variables. Let the edges of 


1 


the plate be pivotally supported. Then, from (21) for k = 0, q = 8.193, it follows p, = MIN Py =P, =P)» = 6.580. 


Two forms of equilibrium correspond to this critical value. Table 1 shows the values of the critical load p’ loss of 
stability of an orthotropic plate with pinched edges compressed along the x axis for different values of the modulus of 
subgrade reaction k and different values of the compressive load along the y axis. Table | shows that with an increase in 
the value of the modulus of subgrade reaction, the critical load of the loss of stability of the compressed orthotropic 
plate also increases. 

Table | 


Dependence of the critical load p’ on the base parameter k at different g 


z 
k q=0 q= 10 q=20 
0 47.455 38.011 5.298 
1,000 112.142 105.841 88.1469 
| 2,000 | 146.99 | 144.252 | 122.039 | 


3,000 164.809 163.198 141.578 
4,000 177.741 176.455 175.106 
5,000 189.002 187.843 186.65 
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Example 2. Let us consider the same orthotropic plate from Example | and direct the x axis across the fibers of the 


face. Then Young's modules along the x and y axes are E, =—-10%kg fom’ E, =1.4-10°kg/cm” respectively; 


Poisson's coefficients — v,=0.46/12 and v,=0.46; shear modulus for the main directions of elasticity 
5 2 er : : ; D,_E, 

G=1.2-10°kg/cm’ . Then, for the problem in dimensionless variables (12), we obtain 6 1; c, D E 12; 
1 1 


D 
C3 = a. = 2.481. Let the boundary conditions (13) be fulfilled — the edges of the plate are freely pinched. As a result 


1 
of numerical calculations using a computer program for solving the problem (12), (13) for this orthotropic plate, 
dependencies are constructed between the load g on the y axis and the critical value p, of the load p on the x axis at 
fixed values of the remaining parameters of this problem. Table 2 shows the values (q, p *) atk=0. 
Table 2 


Dependence between g and p* atk=0 


No. 1 2 3 4 ) 6 
q 0 100 200 300 400 500 
p* 254.388 | 241.867 | 227.789 | 208.778 | 166.462 77.105 


It should be noted that all six points (q, p*) correspond to the form of equilibrium, gradually transforming from the 
form shown in Fig. | a to the form in Fig. 1 c. Figure 1 illustrates the first, third and sixth eigenfunctions corresponding 


to the critical points (q,p*) of Table 1. These eigenfunctions show the forms of equilibrium of the plate after loss of 


stability (forms of loss of stability). 


Fig. 1. Forms of plate stability loss: a — at q = 0, p = 254.388; 
b — at q = 200, p = 227.789; c — at g = 400, p = 166.462 


As a result of numerical calculations using a computer program for solving the boundary value problem (16), (17) 


for this orthotropic plate, parameter values are found k= 200, q=0 for which two eigenfunctions w, and w, 
correspond to the eigenvalue p* = 255.6. This case corresponds to the loss of stability of the plate in two proper forms. 


Figure 2 shows graphs of the forms of stability loss corresponding to this case. 
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Fig. 2. Two forms of stability loss: a— w,; b— w, 


Example 3. Let us consider an orthotropic plate with parameters E, =1.4-10°kg/cm’, E, = 10" kg/cm’; 


Poisson coefficients — v,=0.46 and v,=0.46/12; shear modulus for the main directions of elasticity 


G =1.2-10°kg/cm’. Let the length of the plate be a=5, the width — b=1, the edges x = + — are pinched, and 


1 
the edges y= “2 are free, the parameter — 6 = 5 (the plate can be considered long), modulus of subgrade reaction — 


k =350. Numerical calculation of the solution of the problem (16) with boundary conditions (17) gives the value of the 


critical load p =37.42, the graph of the corresponding eigenfunction is shown in Fig. 3a. Numerical calculation of the 
solution to the problem (23), (24) at v, =0.46 and k=350 for a beam of length a=5 gives p, =39.01, while the 
eigenfunction has the form shown in Fig. 3b. As you can see, the critical load of the plate stability loss coincides with 


the critical load of the beam with an error |p—p,|=|37.42—39.01|~ 0.04. The eigenfunctions of problems (16), (17) 


and (23), (24) have the same number of half-waves. 


b) 


Fig. 3. Forms of stability loss: a— w (x y) , compressed orthotropic long plate; b — wy (x) , compressed beam 


Discussion and Conclusion. The problem of critical loads of an orthotropic plate compressed in two directions 
lying on an elastic base is investigated. The problem is reduced to solving a parametric linear boundary value 
eigenvalue problem. In the case of boundary conditions of a movable hinge support, exact formulas of eigenvalues and 
eigenfunctions are given, and in the case of free pinching of the edges, a variational method is applied in combination 
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with a finite difference method and a computer program for solving the problem is constructed. The following 
conclusions are obtained: 

1. When the compressive load component gq increases along the y axis, the critical value p* of the load compressing 
the plate along the x axis decreases. 

2. If an orthotropic plate is compressed by a load along a direction that corresponds to greater bending stiffness, then 
the critical value of the loss of stability is greater than the critical value of the compressive load acting along the 
direction of lesser bending stiffness. 

3. In the presence of an elastic base, a greater value of the modulus of subgrade reaction k > 0 corresponds to a 
greater value of the critical value of the compressive load. The presence of an elastic base increases the bearing capacity 
of the compressed plate. 

4. In the case of a long orthotropic plate on an elastic base, the critical load of a compressed plate can be calculated 
approximately as the critical load of a compressed beam on an elastic base with an elastic modulus equal to the elastic 
modulus of the plate in the longitudinal direction. 
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